Skyrmion-antiSkyrmion pairs in Ultracold Atomic Gases 
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We study theoretically the dynamics of two-component Bose-Einstein condensates in two dimen- 
sions, which admit topological excitations related to the Skyrmions of nuclear physics. We show 
that there exists a branch of uniformly propagating solitary waves characterised by a conserved 
momentum. These excitations exhibit a cross-over from spatially extended spin-wave states at low 
momentum to a localised "spin- wave droplet" at intermediate momentum; at still higher momen- 
tum, the conhguration evolves continuously into a Skyrmion-antiSkyrmion pair. We discuss how 
these solitary waves can be generated and studied in experiment. 
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Topological solitons have long been a very rich and in- 
teresting subject in theoretical physics, with applications 
in many disciplines. The Skyrmion, first introduced in 
the 1960s as a topological soliton describing nuclei^lj, has 
since been found to be important in condensed matter 
systems: it describes the charged excitations of certain 
quantum Hall ferromagnets^! and lattices of Skyrmions 
have been observed in chiral magnets [J 13] ■ Recent ex- 
perimental advances have allowed the creation and study 
of Skyrmions in multicomponent Bose-Einstein conden- 
sates (BECs)[3. These systems provide a rich area for 
the study of Skyrmions in three dimensions (3D)[1, 0| 
and two dimensions (2D) d, d, [§| , offering the possibility 
of exploring the statics and dynamics of these topological 
solitons in detail. 

A multicomponent BEC consists of two or more 
bosonic species (different atoms, isotopes or hyperfine 
states) . With this additional degree of freedom, the order 
parameter can be expressed as a "spinor" and defines a 
local spin vector field, £^ as a function of position r[5l,[lo|. 
For two spatial dimensions, a Skyrmion is a topological 
soliton in i{r) characterised by a non-zero topological in- 
dex of the map 5^ ~> 5^ (in 3D, the map is 5^ ^ S'^Q). 
A non-zero index arises when the local spin vector field 
rotates by tt from the centre to the outer edge of the soli- 
ton, in such a way that the 2D configuration sweeps over 
the unit sphere. In multicomponent BECs, Skyrmions 
are a type of coreless vortex, analogous to the Anderson- 
Toulouse vortex in ^He-Ajlli.ll^. In a recent experiment, 
2D Skyrmions were created and detected in a spin-2 BEC 
of ^^Rb, using two spatially modulated Raman beams to 
imprint a Skyrmion profile over three hyperfine states [sj. 

In this paper we study the dynamical properties of 
topological excitations in 2D spinor BECs. We focus 
on a two-component BEC, which also admits configu- 
rations with the topology of Skyrmions The dynam- 
ics of a two-component BEC are constrained by both 
the topological index and by a conserved momentum, 
which itself can be defined in terms of the top olog ical 
density [isl. 14 1. Using these conservation laws 15 1 we 



these solitary waves are Skyrmion-antiSkyrmion pairs, 
moving with a velocity perpendicular to the line separat- 
ing their centres, while at intermediate momentum, we 
find a "droplet" of localised spin waves. Recent advances 
in the control and im aging o f BECs with high temporal 
and spatial resolution 1^ 17 1 will permit these dynamical 
topological excitations to be studied in experiment. 

We study the Gross-Pitaevskii energy functional in 2D 
for a two-component condensate of atoms: 
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(Summation convention is assumed throughout, with i, j 
running over the two components and a, /3 over the two 
spatial dimensions.) We neglect the trapping potential, 
and study the properties in the central region of the cloud 
where the density is uniform. Since Skyrmions can be 
imprinted by exciting atoms between different hyperfine 
states, we take mi — mj for simplicity. The interaction 
parameters Uij for a quasi-2D harmonically confined gas 
can be expressed as [18i|: 
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where the s-wave scattering lengths and is the 

oscillator length along the kinematically frozen axis. We 
are guided by the experimentally relevant system of *^Rb 
in the |1, — 1) and |2, 1) hyperfine states [l^, hereafter de- 
noted component 1 and 2 respectively. For these states, 
Uii ~ Ui2 ^ U22 and so stationary solutions vary on 
lengthscales much larger than the healing length, allow- 
ing us to neglect variations in the total density 7]. 

The two-component condensate wavefunction is a 
spinor which we choose to parameterise as: 
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identiiy a branch of uniformly propagating solitary waves 
of the 2D two-component BEC. At large momentum. 
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where po is the total density and 9 E [0,tt], cf) E [0, 27r) 
and e G [0, 27r) are functions of position r. We can then 
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define the local spin as: 



£ ~ X (^X ~ sin (cos 4>x + sin (py) + cos 9z (4) 

where a is the vector of Pauli matrices. Following [lo| . 
we use this to recast Eq. ([T]): 
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where Vg — (2Ve — cos 6'V(j))h/ 2m is the superfluid ve- 
locity and 



Co 



Cl 



C2 



Uii + U22 + 2C/i2 



Un 



4 

f/22 



Uii + [/22 ^ 2[/i2 



(6) 
(7) 
(8) 



The first term of Eql5] corresponds to the non-linear a 
model, previously studied for a ferromagnetfl^. The 
second term is the hydrodynamic kinetic energy of the 
superfluid flow. The cq term is a uniform energy shift, 
while, in the analogy to the ferromagnet, the ci and 
C2 terms can be interpreted as a magnetic field and an 
anisotropy respectively [lot. ci just shifts the chemical 
potentials for A^i and N2, so, like cq, it has no physical 
consequences and we henceforth define the energy (Eq[S]) 

with Cq — Cl = 0. 

The condensate dynamics are described by the two 
coupled time-dependent Gross-Pitaevskii equations: 
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These equations conserve particle number in each com- 
ponent, Ni, as well as the energy. We will look at con- 
figurations where x ~^ Xoo = (1,0) at spatial infinity, so 
that the particle number in component 2 
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is finite, and related to the excitation size. (Since po is 
constant, iVi is automatically conserved when A'^2 is con- 
served.) Two additional important conserved quantities 
are the topological index, Q, and the linear momentum 
Pa-. 
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where o(r) = ^Sap'^aX*'^ PXi is the topological 
density 13| . (The definition of momentum is re- 
lated to the hydrodynamical impulse of classical fluid 
dynamics [20}.) Under the condition x ~^ Xoo at spatial 
infinity, the topological index, Q, is an integer. 



For given Pa, Q and N2, we find the wavefunction 
configurations which minimise the energy E. The time 
evolution of these configurations follows from Eql9]as: 
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where Va and to are Lagrange parameters, introduced to 
enforce the constraints on Pa and N2. It is straightfor- 
ward to demonstrate that this is consistent with trav- 
elling wave configurations, which uniformly propagate 
through the system with a constant velocity Va, while 
the local spin precesses around £z at angular frequency 
UJ. The "energy dispersion" EQ{Pa, N2) is the minimal 
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N and Q, and so: 
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Following Ref. 15, solutions with non-zero velocity ex- 
ist only found for Q^O. (Solutions with Q ^ in- 
clude static Skyrmions, or mutually-rotating Skyrmion- 
Skyrmion pairs.) Henceforth we restrict our analysis to 
the topological subspace Q=0, and drop the subscript 
Q. We shall find that some of the extremal configura- 
tions are localised, and these are therefore the uniformly 
propagating solitary waves 15|, |21 |. 
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FIG. 1. The energy dispersion of a two-component condensate 
for 77 =0.015 and 0.030 and £2=0. The energy for a system of 
free spin waves, Ea llSI is shown for comparison in the limit 
of an infinite system, rj 0. 

We have investigated the discretised energy functional 
E* {Pa, N2) numerically over a square lattice of length 
L. Here, we present results for a system of size 115 by 
115, with the boundaries set to x = Xoo- The scale- 
invariance of the theory requires that under a scaling 
A, E*{P,N2,L) = E*{P\,N2\^,L\)\M- This implies 
that the energy can only depend on the "scaled mo- 
mentum" p = P/ \J N2Poh^ , the "boundary parameter" 
77 = N2/{poL'^) and the "anisotropy" £2 = N2C2m/h^. 
Hereafter we use these dimensionless variables, and ex- 
press the energy in units of h^po/m. The results for an 
infinite system are obtained by taking the limit 77 — > 0. 
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FIG. 2. The local spin vector £{r) projected onto its x — y 
components, and contours of the typical particle density in 
component 2, for p=4.0 {i.e. p > p*) and 77=0.015. For 
clarity, only one spin vector in five is plotted. The direction of 
the velocity and momentum of the solitary wave is indicated. 
Both maps show a localised structure resembling free spins, 
therefore termed a spin-wave droplet. 
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FIG. 3. As in Fig[2] but now for p—8.5. Both maps show a 
localised structure resembling a Skyrmion-antiSkyrmion pair. 



Figure [T] shows the minimum energy as a function of 
p for C2—0, at two values of the boundary parameter, rj. 
Although this behaviour is still weakly dependent on 77, it 
is clear that E* is tending to a smooth function as 77 — > 0. 

We have analysed the extremal configurations and find 
that as 77 — >■ 0, there is a transition from spatially ex- 
tended configurations at small p to localised structures 
at large p. A study of the polar angle at the centre of 
the system, 6c, shows that for p > 2.4, 9c tends to a 
finite value in the limit of an infinite system indicating 
localised solutions. Below p ~ 2.4, the system contains 
delocalised configurations with 9c extrapolating to zero 
as 77 — > 0. From this, we identify a transition at p*=2.4, 
consistent with that found for the ferromagnetjlSj. The 
energy of delocalised states below p* can be understood 
within a linearized continuum theory [l5j|. which gives: 



Esw = y + Ti" 77 
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corresponding to a system of length L filled with N free, 
non-interacting spin waves. Numerically, we find that E* 
approaches this form as 77 — >■ 0, and the particle densities 
in each component match those for free spin waves. 

Above p* , the spatially localised configurations form 
the branch of uniformly propagating solitary waves. Typ- 
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FIG. 4. The energy dispersion of a two-component condensate 
for C2 = —2.0, 0.0 and 2.0 for 77 =0.015. Energy is measured 
relative to the uniform system x ~ Xoo- The value £2 = —2.0 
is below the critical value £2 = —1.5, and so solitary waves 
exist at all momenta, p, for 77 -> 0. Inset: Enlargement of 
the bifurcation region for £2 = 2.0, within which we find two 
branches of excitations for the same momenta. 



ical spin configurations and particle densities in compo- 
nent 2 are shown in Figsl2]&[3]for intermediate and high 
p respectively. At intermediate momentum, the structure 
resembles free spin- waves with a finite spatial extent and 
so is referred to as a "spin- wave droplet" . 

At high momentum, the local spin configuration 
matches that for a Skyrmion-antiSkyrmion pair. From 
Eq ll2l the scaled momentum, p, is related to the sep- 
aration of a Skyrmion-antiSkyrmion pair, r, by p = 
2rTTy^ P0/N2. For wide separations {i.e. large momenta), 
the dominant contribution to the energy, E*, is from 
the kinetic energy of the superfluid flow. The prop- 
erties of this regime therefore differ from those of the 
ferromagnet previously studied [15J. Outside the core 
regions of size vq (wherein £ describes a Skyrmion or 
anti-Skyrmion), the superfluid flow is that of a vortex- 
antivortex pairj22|. The scaled energy, E* , then includes 
the term 27rlnfr/rn)(20|. The energy dispersion curve 
at high p approaches this result continuously as 77 — 
and the velocity, v, of the Skyrmion-antiSkyrmion pair 
(EqfT4|) tends to v — Ti/mr. 

We have investigated how the anisotropy, C2, affects 
these results. Figure |4] shows the minimum energy as a 
function of p for £2 = —2, and 2, calculated with respect 
to 

X — Xoo- Negative values of £2 represent an "easy- 
axis" anisotropy, preferring spatial separation of the com- 
ponents. The ground state is one in which all atoms are 
in a single component, and the extremal configurations 
correspond to those discussed above. As the anisotropy is 
decreased below zero, the critical scaled momentum, p*, 
shifts down, until below a critical value, C2, the solitary 
waves exist for all values of momentum. Applying the ra- 

4c2 from 



dial ansatz of Ref. |15|, we find that p* = -v/5-8 - 
the solution of a ID nonlinear equation [23|. Hence, the 
transition point, p* , vanishes below a critical anisotropy, 
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C2 = —1.5, which is consistent with our numerical results. 

A positive £2 represents an "easy-plane" anisotropy, 
preferring component mixing. The ground state will be 
an equal mixture of the two components and, instead 
of Skyrmions, the coreless vortices are half-Skyrmions or 
merons (also called Mermin-Ho vortices) . where the local 
spin rotates by 7r/2 radially outwards [Sll 1 2ll24l| ■ 

However, even for positive £2, Skyrmions can be cre- 
ated and can be dynamically stable: starting from a sin- 
gle component BEC, such as x — Xoo (with N2 — 0), ex- 
citation into a state with non-zero N2 and P can lead to a 
solitary wave which cannot relax to the ground state due 
to the conservation laws. The energy dispersion curve 
for these solitary waves is shown in Fig HI We find that 
a positive value of £2 increases the critical scaled mo- 
mentum, p* . Furthermore, it introduces a bifurcation 
region, where two excitation branches exist for the same 
momenta. From our numerical results, we find that exci- 
tations in the higher momentum branch always have the 
form of the Skyrmion-antiSkyrmion pair described above. 

Our results show that there is an interesting cross- 
over from spatially extended spin-wave states at small 
momentum to localised spin-wave droplets at inter- 
mediate momentum. As momentum is further in- 
creased, the spin- wave droplet evolves continuously into 
a Skyrmion-antiSkyrmion pair, forming a branch of uni- 
formly propagat ing solitary waves. Recent experimen- 
tal advances [5|, |l6| should allow the imprinting of a 
Skyrmion-antiSkyrmion pair on a quasi-2D multicompo- 
nent BEC and the study of the dynamical evolution with 
high temporal and spatial resolution. Varying the separa- 
tion of the imprinted pair varies the momentum, allowing 
the dynamics of topological solitons to be investigated 
experimentally for the first time. 

The anisotropy, £2, controls the position of the cross- 
over and the type of extremal configuration. For the 
1) and |2, 1) states of ^''Rb, the scattering lengths 
are: 022 = 95.00ao, an = 100.40ao and 012 = 97.66ao 
where qq is the Bohr radius [2^, giving a small positive 
value of C2 (Eq|51). For a density of tt-sd — lO^^cm"'^, 
the anisotropy healing length is ^2 — 30/im. Anisotropy 
£2 is important only for excitations with C ^ ^2, where 
£ is the scale over which the density of component 2 
varies. Recent experiments have demonstrated a spatial 
resolution on scales as small as 0.5/im[ll], and so would 
allow an exploration of regimes of both weak and strong 
£2 (for £ < 30/im and C > 30/Km respectively) [26] . 

The solitary waves travel through the system with a 
uniform velocity, v, while the local spin precesses around 
£z with angular frequency, uj. For an excitation with size 
C ~ 3/im at p = 8 for ^'^Rb, our results typically show 
that V ~ 0.2mm s~^ and w ~ 0.3rad ms~^ (EqlTi]). Thus, 
the temporal resolution required to resolve w is ~ 20ms 
and it will take ~ 50ms to move a distance of 10/xm, both 
of which are within current experimental capabilities [T?]. 

The creation and characterisation of the solitary 



waves predicted here will allow subsequent experimental 
studies of their scattering properties. It would be espe- 
cially interesting to investigate the general prediction 
coming from high energy physics, that pairs of (related) 
2D solitary waves of equal and opposite momenta should 
scatter at right angles after coUiding head-on(27i l28j. 
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